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UNIVERSITY QI ILLINOIS 
"DIGITAL COMPUTER 

LIBRARY ROUTIN E SA 3 - 155 
Natural Logarithm (DOI or SADQl) 

Subroutine (entered -with 8J order from Library Routine Al 
15 milliseconds 
30 

The subroutine determines the natural logarithm of an 
arbitrary (positive) number 'which lies in the floating 
accumulator in the standard form Hue. (l/lO < A < l) 
where A is the fractional part]. Control is returned 
after computation to the next floating order. 
At all times the accuracy of the floating point routine 
is available (provided the argument "was in its standard 
for initially). For arbitrary exponent there is a 
maximum error of + 1 in the 9"th decimal place. 
0, 1, 2 

Enter vlth an 8J order from a set of interpretive orders 
with the argument In standard form in the floating accumu- 
lator. After calculation, control returns to word 19 of 
Library Routine Al. 

The result appears in the floating accumulator in 
standard form» 
S3' 
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Floating accumulator 



J 



S k address of first word of Library Routine Al 

If the argument is a negative number or the program will 
•stop on a division hang -up. 

RT: 10/13/60 



DAT E November 9, 195^ 
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"ft 5 ? — 

L0 4184 
40 2F 

5^ 20L. 
75 2F 
S5 F 
40 F 

F5 130S4 
42 1S3 
L5 S3 
L4 19L 
kO IF 
L5 S3 
FO 19L 
66 it 
S5 22L 
40 IF 
7 J IF 
kO 2? 

L5 7L 
42 11L 
50 21L 

7J 2F 
L4 ( ) L 
40 S3 
50 S3 
F5 11L 
40 I1L 
LO 18L 
36 11L 
7J IF 
LO SOL 
L4 F 
L4 f 
ho S3 
26 19S4 J 
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yram- p/400 In 10 



Adjust exponent for scaling 
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Cdsputs 1/200 la. y 



Cospute a/200 (k y » In flo'+ p In 10 j 
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LT 2F 
L4 30L 




End constant 
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00 F 








oo 15811 38830 OOJ 




1/2 yio 
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00 F 

00 57564 62732 J 




1/200 logJlD 
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00 F 

00 32085 28548 J 




c 1 /200 = 0.0052085285^8 
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NO F 

00 1*9832 89053 80J 




C 15 /200 = -0,001671094620 
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00 F 

00 17477 06534 J 




c /200 = 0.001747706534 
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00 F 

00 68285 8823J 




c. , /200 = .000682858823 1 
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00 F 

00 Il4l0 I6968 J 




0^200 = 0.001141016968 
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00 F 

00 14263 85504 J 




C ? /200 = 0.001426385504 
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OQ F 
iOO 20000 79233 J 




C.-/200 = 0.002000079233 
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00 F 

00 33333 32235 J 




CJ200 = 0.003333332235 
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00 F 
-00 10000 00000 OJ 




C.j/200 = 0.010000000000 
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